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Abstract 



Wc consider a stochastic fluid queue served by a constant rate server and driven 
by a process which is the local time of a certain Markov process. Such a stochastic 
system can be used as a model in a priority service system, especially when the time 
scales involved are fast. The input (local time) in our model is always singular with 
respect to the Lebesgue measure which in many applications is "close" to reality. We 
first discuss how to rigorously construct the (necessarily) unique stationary version of 
the system under some natural stability conditions. We then consider the distribution 
of performance steady-state characteristics, namely, the buffer content, the idle period 
and the busy period. These derivations are much based on the fact that the inverse of 
the local time of a Markov process is a Levy process (a subordinator) hence making the 
theory of Levy processes applicable. Another important ingredient in our approach is 
the Palm calculus coming from the point process point of view. 

Keywords: Local time, fluid queue. Levy process, Skorokhod reflection, performance 
analysis. Palm calculus, inspection paradox. 

AMS Classiflcation: 60G10, 60G50, 60G51, 90B15. 

1 Introduction 

This paper extends the results of Mannersalo et al. [I3] who introduced a fluid queue (or 
storage process) driven by the local time at zero of a reflected Brownian motion and served 
by a deterministic server with constant rate. The motivation provided in [13] is that the 
system provides a macroscopic view of a priority queue with two priority classes. Indeed, 
in such a system, the highest priority class (class 1) goes through as if the lowest one 
does not exist, whereas the lowest priority class (class 2) gets served whenever no item of 
the highest priority is present. In telecommunications terminology, class 2 only receives 
whatever bandwidth remains after class 1 served. As argued in |13j . if the highest priority 
queue is, macroscopically, approximated by a reflected Brownian motion, the lowest priority 
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queue is driven by the cumulative idle time of the first one, which is approximated by the 
local time of the reflected Brownian motion at 0. 

In view of Internet networking applications, such as the service provision amongst several 
classes of service (e.g. streaming video and expedited data), the fluid or macroscopic model 
is thus quite appropriate for obtaining a better picture of the situation and for performance 
analysis and design. From a mathematical point of view, the model is a rare example of a 
non-trivial fluid queue whose performance characteristics (such as steady-distribution) can 
be computed explicitly. If, in addition, we take into account the heavy-tailed nature of 
traffic on the Internet, it seems reasonable to consider a Levy process as a model for class 
1 queue. This provides motivation for studying a queue whose input is the local time of a 
reflected Levy process. 

More generally, let X be a Markov process and L its local time at a specific point. The 
fluid queue driven by X refers to the stochastic system defined by 

Qt = QG + Lt-t + Iu t>0, 
where > for all t > 0, and / is a non-decreasing process, starting from 0, such that 

/ 1(Q, > Q)dls = 0. 

^0 

Thus, Q is obtained by Skorokhod reflection and / is necessarily given by 

It = - inf [{Qq + Ls- s) 

0<s<t 

see [H]. By considering, instead of 0, an arbitrary initial time, we can define a proper 
stochastic dynamical system (see Appendix A for details) which, under natural conditions, 
admits a unique stationary version. To this we refer frequently throughout the paper. 

We remark also that in Kozlova and Salminen [TO] the situation in which X is a general 
one-dimensional diffusion is analysed. Moreover, Sirvio (nee Kozlova) [16] studies the case 
where L is constructed as the inverse of a general subordinator (without specialising the 
underlying process X). 

This paper follows ideas which were developed in [16] in the context of refiection of 
the inverse of a subordinator. However, (1) it connects the abstract framework with the 
case where the subordinator is the local time of a reflected Levy process (motivated by 
applications in priority processing systems) (2) it uses, as much as possible, a framework 
based on Palm probabilities (see, in particular. Section 4-Theorem [3] and Section 5-Lemma 
[2]), and (3) explicitly discusses the possible types of sample path behaviour of the process 
of interest (Q); for illustration, see Figures 1, 2 and 3: Figure 1 concerns the case where Q 
has continuous sample paths but with parts which are singular with respect to the Lebesgue 
measure. Figure 2 concerns the case where Q has paths with isolated discontinuities (positive 
jumps) and linear decrease between them. Figure 3 concerns the case where Q has absolutely 
continuous paths. The cases are exhaustive. Despite the wide variety of sample paths 
(depending on the type of underlying Levy process Y), the mathematical framework and 
formulae derived have a uniform appearance. 

The paper is organised as follows. In Section [2] we construct the stationary version of 
the underlying (background) Markov process X. In Section [3] we construct the stationary 
version of stochastic fluid queue with input the local time of X, based on the stationary 
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version of X. In Section H] we derive the stationary distribution of the buffer content and 
present a number of examples. In Section [5] we examine the idle and busy periods, and, in 
particular, characterise the distributions of their starting and ending times. The analysis 
is carried out first under the condition that the time at which the system is observed is a 
typical point of and idle or a busy period. Finally, the distributions of typical idle and busy 
periods are derived. 



2 The background Markov process and its local time 

We first construct the underlying Markov process X which models the highest priority 
class. This process will be taken to be the stationary reflection of a spectrally one-sided 
Levy process 

Y = {Yt, t G M) 

with two-sided time and Yq = (see Appendix B). A Levy process is called spectrally 
negative if its Levy measure 11 satisfies 

n((-oo,0))>0 and n((0, +cx))) = 0, 

and spectrally positive if 

n((-oo, 0)) = and n((0, +oo)) > 0. 

Clearly, if Y is spectrally positive then —Y is spectrally negative, and vice versa. To avoid 
trivialities we shall throughout assume that 

\Y\ does not have monotone paths 

which rules out the cases that Y is an increasing or decreasing subordinator. 

We also discuss the characteristics of its local time at 0. Appendix A summarises the 
notation and results on the Skorokhod reflection problem and its stationary solution. Ap- 
pendix B summarises some facts on Levy processes with one sided jumps indexed by M. We 
will throughout denote by P a probability measure which is invariant under time shifts, and 
by Px the conditional probability measure when Xq = x. We define the Laplace exponent 
of a spectrally one-sided Levy process as a function ipy ■ M+ — > M given for > by 



log Ee^^^'+^-^'\ if Y is spectrally negative 
log -Ee~^^^*+^~^*\ if Y is spectrally positive. 



Thus we insist that ipy be defined on M+j and define its right inverse 

<^Yiq) := sup{e > : Vy(^) = q}, q> 0. (1) 
We use also the notation 

Yt := sup Ys and Y^ := inf Ys, t > 0, 

0<s<t 0<s<t 

and recall the duality lemma for Levy processes (see, e.g., Bertoin [21 p. 45]): 

{Yt - Y^t_,y : 0<s<t} = {Ys : 0<s<t}, (2) 
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which is equivalent with 



where = means equahty in distribution. Hence, 

sup (yt-y(t_,)_) =Ff, 

0<s<t 

Yt-Y,^Yt. (3) 

In this paper, we wih mainly study the Levy process Y with the time parameter taking 
values in the whole of M (see Appendix B). The Skorokhod reflection mapping associated 
with {Yt : t € M} is defined (see Lemma [8] in Appendix A) via 

Xt:=ntY:= sup (Yt-Y,), t € M. 

— oo<s<t 

In the remaining of this section we give conditions for the existence of the stationary proces^ 
TZY = {TZfYji G M), compute its marginal distribution, and define the local time process L 
of TZY at which will be used for the construction of the fluid queue. A few words about 
the definition of L are in order. We adopt the point of view the L is a stationary random 
measure on (M, <^), i.e. 

L{s,s + t]= L{0,t]oes, t>0, sGM, 

where is the shift on the canonical space (see Appendix B), which regenerates together 
with X at each point t at which Xt = 0. It is known that L is a.s. continuous if and only 
if the point x = is regular for the closed interval (— cxo,0] for the process X, and this is 
equivalent to inf{t > : ^ < 0} = 0, Po-a.s. Furthermore, L is a.s. absolutely continuous if 
and only if, in addition to the above, the point x = is irregular for the open interval (0, oo) 
for the process X, and this is equivalent to inf{t > : If > 0} > 0, Po-a.s. If L is a.s. 
continuous then it is not difficult to attach a physical meaning to it as a cumulative input 
process to a secondary queue. For mathematical completeness, we shall also consider the 
case where L is a.s. discontinuous, in which case it can be shown to have a discrete support. 
In the continuous case, L can be defined uniquely module a multiplicative constant. We 
shall make the normalization precise later. In the discontinuous case, there is more freedom; 
however, insisting that its inverse be a subordinator, we are left with only one choice. The 
discontinuous case appears only once below and the construction of L is discussed there. In 
all cases, the support of the measure L is the closure of the set {t G M : = 0}. 

Associated to the measure {L{B),B G ^) we can define a cumulative local time process, 
denoted (abusing notation), by the same letter, and given by: 



The right-continuous inverse process is 



L{0,t], t>0, 
-L(t,0], t<0. 



_^ (mi{t > : L(0, t]>x}, x>0 
' |sup{t<0: L{t,0]<x}, x<0. 

In case L is P-a.s. continuous, the process {L~^,x G M) has, under Pq, independent in- 
crements and a.s. increasing paths (i.e. it is a, possibly killed, subordinator). This is an 
additional requirement that needs to be imposed when L is not P-a.s. continuous. 
Throughout the paper, we denote by Px the probability P, conditional on Xq = x. 



^That this process is Markov is easy to see due to the independence of the increments of Y . 
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2.1 Stationary reflection of a spectrafly negative Levy process 

Suppose that the process 1" is a spectrally negative Levy process with non-monotone paths; 
see expression ([53]) . 

Proposition 1. Let Y = {Yj : t € M} he a spectrally negative Levy process with two-sided 
time. Assume that its Laplace exponent ipviG) = logEe^^^ , 6 > 0, is such that i/'y(0-|-) < 0. 
Then 

X = {Xt := TZtY : t G R} 

is the unique stationary solution of SDS ( the Skorokhod dynamical system, see Appendix A ) 
driven by Y. The marginal distribution of X is exponential with mean l/<J>y(0). 

Proof. Since E[Yt^i — Yt] = V'y (0+) < 0, existence and uniqueness of the stationary solution 
is guaranteed by Corollary [2] of Appendix[Al That <&y(0) > is a direct consequence of the 
definition of <I>y (see ([T])). To derive the marginal distribution of X consider for (3 > 

Since Y is spectrally negative its over all supremum, supg>oys, is exponentially distributed 
with mean l/<l>y(0) (see, e.g., Bertoin |2j p. 190] p. 190 or Kyprianou [HI p. 85]). □ 

In view of Proposition [H we assume that 

ij'Y{0+) G [-00,0), 

which is equivalent to 

$y(0) > 0. 

It is easily seen that, for a non-monotone spectrally negative Y, a necessary and sufficient 
condition for continuity of L is that Y has unbounded variation paths. This is further 
equivalent to: a > or J^-^ \y\Il{dy) = oo. 

In the alternative case, when the paths of Y are of bounded variation, the number of 
visits of Y to its running infimum forms a discrete set. So n{s,t] := YIugr ■^(■^'u. — 0) 
finite for all — oo < s < t < oo. Let nt := n(0, t] if t > 0, and nt := —n(t, 0] if t < 0, and let 
{tj,j G Z) be a collection of i.i.d. exponentials with mean 1, independent of Y. We adopt 
the following construction for L. 

-^(•S)^] = — OO < s < t < oo. 

na<i<nt 

In both cases, the process {L~^,x G M) is a subordinator under Pq, with Lq ^ = 0, Pq-b-.s. 
If L is continuous, this property is immediate from the definition of L. If L is discontinuous, 
L^^ has independent increments due to our choice of the exponential jumps of L. Since 
"0^^(0-1-) < 0, we have Yt Too, as t ^ ±00, P-a.s., and this implies that Lt ±00, as 
t — > ±00, P-a.s. Thus, it is not possible for L~^ to explode for finite x. 

Regardless of the continuity of the paths of L, we always have the following: 

Proposition 2. Let Y and X be as in Proposition [I] and L the local time at of X. Then 
the local time L can be normalized to satisfy 

^o[e-'?-^^'] =e-^''/*^(''), q>0, (5) 
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and, moreover, 



EqL-^ = x/^y{0), xeK, ELt = t^Y{0), t€R. (6) 

Proof. For the Laplace exponent of in ([5]), we refer to Bingham [3l p. 731] (a result 
due to Fristedt [5]) and Kyprianou [11] and Kyprianou and Palmowski [12]: The "ladder 
process" {L~^ , X^-i),x > 0) is a Levy process with values in and Laplace exponent 

log Eo[e . ] = ^^(,)_^ 

obtained by Wiener-Hopf factorisation for a spectrally negative process; see Bertoin [21 p. 
191, Thm. 4]. Setting /? = we obtain £'o[e~"^i ] = e~"/*^'("), as claimed. From this, we 
obtain EqL~^ = x/^'y(0), by differentiation. Using the strong law of large numbers, we 
have liuix^oo / X = l/^>y(0), Po-a.s., and, given that > 0) is the right-continuous 

inverse function of {Lt,t > 0)-see (jH, we have limt-,oo Lt/t = <l>y(0), Po-a.s., and P-a.s. 
Since L is a stationary random measure, we have ELt = C t for some constant C. Hence, 
we immediately have that C = <I>y(0), and this proves ([6]). □ 

We shall later need the P-distribution of the random variable 

D := inf{t >0: Xt = 0}. (7) 

Since Xq is exponential with rate <^y(0), we have ^(-'^o > 0) = 1. So if < t < D, we have 
Xt = Xq + Yt, P-a.s. Therefore 



Let 



D = M{t > : Yt< -Xo}, P - a.s. 
r_a; := inf{t > : Yt < -x}. 



Therefore, 

/•oo 

E[e-^^] = / P{Xo e dx) E[e-^^- 
Jo 



'<I>y(0)e-*^W^ (^Z^'^){x)-^^W^'\x)^ dx 
V'y(^y(O))$y(0) -0$y(O) 



'^'^^^^ -^Yie) ci>y(o) [i;Yi^Ym)-e] 



$y(0)' 

where the overshoot formula (j6ip given in the Appendix B (see also Bingham [3l p. 732]) 
and the Laplace transforms (I56p . (I57p . for the scale functions W^^\ Z^^\ were used. 



2.2 Stationary reflection of a spectrafly positive Levy process 

We can repeat the construction in the subsection above for a spectrally positive Levy process 
Y with non-monotone paths. We shall be using the formulae of Appendix B with —Y in 
place of Y. 
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Proposition 3. Let Y = {Yt : t G M} be a spectrally positive Levy process with two sided 
time and Laplace exponent tpy- Assume that its Laplace exponent ipvid) = logEe~^^^ is 
such that iP'y{0+) > 0. Then the process {Xf := TZfY : t G M} is the unique stationary 
solution to the SDS driven by Y. The stationary distribution of X is given for P >0 by 

t—^00 t— ►OO 

Proof. Notice that in this case, by the assumption on ■i/'y, 

E[Yt+i - Yt] = -M0+) < 0, 
and, hence, Y drifts to — oo. Let Z := —Y. Clearly, Z is spectrally negative, 

= Yt — y.t — — -^t) 

and Z drifts to +00. The classical result due to Zolotarev [18] (see also Bingham [3] Propo- 
sition 5 p. 725) says that the stationary distribution of X is as given in Q. □ 

We shall therefore assume that 

V'y(o+) e (0,00). 

Hence $y = tpy^ and so 

^'y{0+) = 1/V'y(0+) € (0,00). 

It is easily seen that, starting from 0, the process Y hits (— oo,0] immediately, P-a.s., 
and this ensures continuity of the local time L. Moreover, we may and do normalize L so 
that 

L{s,t] = - inf Y{s,u]. (10) 

s<u<t 

The continuity of L implies that 

{L-i : X > 0} = {r_, : X > 0}, (11) 

where T-^ := inf{t > : Yf < -x} = inf{t > : Zf > x}. Note that iL~^,x £ M) is 
a subordinator under Pq, with Lq^ = 0, Po-a.s. Furthermore, since Yt drifts to ^00 as 
t — > ±00, L~^ is proper (not killed). 

Let us briefly comment on the special case where, starting from Xq = 0, the interval 
(0, 00) will be first visited by X at an a.s. positive time. It is known [21 Ch. 7] that this 
occurs if and only if Y has bounded variation, i.e. 

Y{s,t]=dY{t-s)+ f [ yr]{du,dy), (12) 

J{s,t] J(0,oo) 

where dy is the drift, and Jq y Il{dy) < 00. Since we exclude the case where Y is monotone, 
we must have dy < 0. In that case, with ()10p as the definition of L, it is known that for all 
s < t, 



L{s,t] = \dy\ [ l{Xu = 0)du. 

J s 



(13) 



A rewording of the first part of Lemma [10] in Appendix B gives the first part of the 
following proposition. 
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Proposition 4. Let Y be as in Proposition and L the local time at of X. Then 



£;o[e-«-^-'] =e-*^(«)^, q>0. (14) 

Moreover, 

EoL'^ = x^y{0+), x>0, ELt = t'il)Y{0+), t>0. (15) 



Proof. Formula (I14p follows from (llip and the well known characterisation of the distribution 
of the first hitting time t^, see, e.g., Bingham [3] p. 720. By differentiating (jl4p . we obtain 
the first part of (jlSp . and using an ergodic argument-as in the proof of the second part of 
([6])-we obtain the second part of (fT5|) . □ 

We now compute the P-distribution of D = inf{t > : Xt = 0}, by arguing as earlier: 
we have D = mi{t > : —Yt > Xq}, and, since —Y is spectrally negative, we use the hitting 
time formula (j60p of Appendix B to obtain 

^[,-», _ ^[,-»,(.).v„, _ _2vl|_ _ 

where we also used ([9]) and the fact that ^y°'4'y is the identity. 



3 Construction of (the stationary version of) the fluid queue 
with local time input 

We wish to construct a fluid queue driven by 

L{s,t] = L{s,t] - {t- s), 

where L is the local time at zero of the Markov process X. The process X is a stationary 
Markov process which is the reflection of a spectrally negative (Section 12. ip or a spectrally 
positive (Section 12. 2p Levy process. In either case, L is a stationary random measure with 
rate (see dSD and ([l5])) 



j^'y(O), if y is spectrally negative 
^ = EL{s, s + 1\ = < . (17) 

I '!/'y(0+)) if Y is spectrally positive 

The fluid queue started from level x at time is, as explained in Appendix [Al the process 

(7^o,^^(x), t > 0). 

From Corollary [2] in Appendix!^ we immediately have: 

Theorem 1. If fj, < 1 there is a unique stationary version of the fluid queue driven by L 
and is given by 

Qt = jZtL= sup (Lt-K), tern. (18) 

— oo<-u<t 

Thus, the following assumptions will be made throughout the paper: 
[Al] n < 1 

[A2] In case Y is spectrally positive and of bounded variation then dy < —1 



8 



Assumption [Al] is so that we can construct a stationary version of Q (as in Theorem [T]) . 
If Y is spectrahy positive and of bounded variation non-monotone paths then its drift dy 
must be negative. If, however, |(iy| < 1 then-see (fT3l) -L(s. t] < t — s for all s < t and so Q 
will be identically equal to 0. 

Physically, we think of Q a stationary fluid queue whose cumulative input between times 
s and t is L(s,t] and whose maximum potential output \st — s. Unlike X, the process Q is 
not Markovian. However, since Q has been built on the probability space supporting X, it 
makes sense to consider, for each x > 0, the probability measure Px defined as P conditional 
on {Xq = x}. 

Since L is a random measure allows us to consider the Palm distribution with respect 
to it, namely the probability measure defined by 



PL{A) = ti-^E[ lAoeMdt). 

J(0.1] 



'{0,1] 

Since Pl is also given as the value of the Radon-Nikodym derivative 

^^ElAL{dt) 



dt 



t=Q 



it follows that P^ expresses conditioning with respect to the event that t = is a point of 
increase of L. But L is the local time of X at zero. Therefore, the support of L is the set 
of zeros of X. We thus have 

Theorem 2. The Palm measure Pl coincides with Pq. 

This observation allows us to use the formulae of Appendix [Al involving Palm probabil- 
ities. 



4 Stationary distribution of the fluid queue 

We are interested in computing P{Qt G •), a probability measure which is the same for all 
t. We will use three properties of L. First, duality, i.e. that Y has the same distribution 
when time is reversed, see ([2]), implies that 

(L(0,t], t > 0) = {L{t,0), t < 0), under P and under Pl. 
Second, the process 

L"^ = inf{t > : L(0, t]> x}, x>0 

is a subordinator under Pq. Third, the Palm measure Pl coincides with Pq. 

Recall that V'y(^) has been defined as log -56^^^'+^"^*^ when Y is spectrally negative and 
as log i?e~^'-^'+^~^'^ when Y is spectrally positive. The reason is that it is customary to 
have ^ > in both cases. Recall also that $y(g) = sup{0 > : iPy{G) = <?}• The stationary 
distribution of Q will be expressed in terms of ipY- For earlier works on this problem, we 
refer to [15] for diffusion local times and [16] for the inverse of a general subordinator. The 
present formulation in Theorem [3] is in particular tailored for the local time of X. The proof 
makes use of the Palm probability which is a new ingredient. 

First, since we allow discontinuous local times, the following simple lemma is needed. 
We omit the proof. 
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Lemma 1. Let M := MU {+00, —00}, ///i : M ^ M is right- continuous and non- decreasing 
then h~^{x) := inf{t : h{t) > x}, x € M, is right- continuous non- decreasing, h^^ : M — > M, 
and, for all t G M, 

h~^{h{t-)-) < h~\h{t)-) < t < h~\h{t-)) < h-\h{t)). (19) 
Furthermore, = h. 

Theorem 3. (i) If X is the reflection of a spectrally negative Levy process Y with i/;y(0+) < 
0, 7pY{l) > then 

Po(Qo >a)= e-'^^(l)^ P(Qo > a) = $y (0)e-^^(l)^ a > 0. 

(a) If X is the reflection of a spectrally positive Levy process Y with < 'i/'y(0+) < 1 and 
dy < — 1 in the case of bounded variation, then 

Po{Qo> a) =6-^*", P(Qo > a) = ^y(0+)e-^*^ a > 0, 

where 6* > is defined by ^y(0*) = 6* . 

Proof. By the construction of Q and duality, we have P{Qo < a) = P(sup„>o(-Z^u — u) < a). 
At this point we note that our assumptions imply that the process {Lt — t:t>0} does not 
have monotone paths. The event {supj>o(-C't — t) < a} can be expressed in terms of L~^: 

{sup(Lf — t) < a} = {sup(x — L~^) < a}. 

t>0 x>0 

To justify this (recall that in case (i) L is not necessarily continuous), we first assume that 
Lt < t + a for all t > 0. Hence L^^-i < L~^ + a for all x > 0. Since (Lemma [1]) 

Lj^-i > X, X > 0. 

we have x < L~^ + a for all x > and thus sup^>q(x — L~^) < a. Assume next that 
X — L~^ < for all x > 0. Then limgjo -^^-e = -^x- > x — a for all x > 0. Therefore, 
L^L^)- >Lt-a for aU t > 0. Since (by Lemma [1] again) 

L-i)_ <t, t>0, 

we obtain L7^^-^_ > t — a for alH > and this gives sup^>o(it — t) < a. We first compute 
the Po-distribution of Qq: 

PoiQo > a) = Po(sup(x — L^-^) > a), a > 0. 

x>0 

Under Pq, the process 

{A^ := X - L"^ : x > 0}, (20) 
is a spectrally negative Levy process with bounded variation paths. Letting 

aa := inf{x : A^; > a}, 

and applying Lemma [TO] of Appendix B, we have 

Po(Qo > a) = PoicJa < 00) = lim^o^e"''^''] = lim e'^'^^'?)'' = e-**^")". (21) 
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The function $A(g) is given by 

$A(g) = sup{e > : ^a(0) =q}, q>0 

(see (j55|) ). where 

V'aW =logii;oK^i] =^ + logSo[e-'^^']. 

If Y is spectrahy negative, then we use Proposition [2] for an expression for log Eq [e~^'^i ] . 
If Y is spectrally positive we use Proposition [H We obtain: 

\ G — ^ ^iQ\ , if Y is spectrally negative, , , 

I — $y(6'), if Y is spectrally positive. 

In both cases, Qq is exponential under Pq. with parameter $a(0), which has different value 
in each case. Let be the rate of L (see pT|) ). Using equation ([5T|) . we have 

^(Qo>a)=/^^'L(Qo>«) = /U^'o(Qo>a)=^e-**(°)^ a > 0. (23) 

The proof will be complete, if we show that 



'^a(O) 



'(/'y(l), if Y is spectrally negative, 
0*, if y is spectrally positive. 



Note that ^a(O) is the positive solution of il^\{9) = 0. If y is spectrally negative, we see, 
from the first of ([22]), that V'a(^) = iff ^yiP) = 1 and, by the definition of $y, the latter 
is true iff 6* = Vy(l)- Thus, $a(0) = V'y(l)- If Y is spectrally positive, 'i/'A(6') = iff 
^y{0) = ^ iff = Vy(^). □ 
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Table showing the basic characteristics of the system in both cases 



Y spectrally negative 



Vy(9) 




<fy(0) 1 



-xq/<S>Y{q) 



fi = rate of L = <I>y(0) 
Po(Qo > a) = e-^*'^ 
r = V'y(l) 

P{Qo = 0) = 1 - $y(0) 




fi = rate of L = ^y(0+) 
Po(Qo > a) = e-'^*" 

P(Qo = 0) = 1 - ^^.(0) 



4.1 Example 1: Fluid queue driven by the local time of a reflected Brow- 
nian motion 

Consider y to be a Brownian motion with drift (see also [15], [13], |10j): 

Yt := aBt - fit, t £ R, 

where a > 0, fj, > 0. Here B = {Bt,t € M) is a standard Brownian motion with two-sided 
time. In other words, {Bt,t > 0), {B-t,t > 0) are independent standard Brownian motions 
with Bq = (although specification of Bq does not affect the results below). The Levy 
measure here is 0. Consider Y as in Section [2.21 and let 



V'y(6l) = log£;[e 



-evy 



1 



+ 0>O. 



Define 



Xt=ntY= sup {a{Bt-Bs)-fi{t-s)), te 

~oo<s<i 



Lemma [3] gives the distribution of Xq under P: 
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i.e. exponential with rate liija^ . Let L be the local time at zero of X. The rate of L-see 
p!7|) -is '(/'y(0+) = \x. Assume /x < 1. Let iv(s,t] = — {t — s) and let Q be defined by 

Qt = 7^^L, t G R. 

Theorem [3] gives the distribution of Qq under P: 

P{Q^ >a)= tP'Y{0+)e-^''' = /ie-2(i-^)«/"', a > 0. 

Here, 9* was found from = 0*. Thus Qq is a mixture of an exponential with rate 

2(1 — /u)/(T^ and the constant which is assumed with probability fi. 

4.2 Example 2: fluid queue driven by the local time of a compound Pois- 
son process with drift 

Suppose that, for a > 0, 

Yt = St- at, t G R, 

where S is a compound Poisson process with only positive jumps, jump rate A and jump 
size distribution F. For simplicity, we take F to be exponential with rate 6 > 0, i.e. 
F{dx) = 6e-^''dx. Then 

,Py{0) = log£;e-^(^'+^-^^) =ae-X [ (1 - e-^'')F{dx) = aO - ^ > 0. 

J[0,oo) + 

The assumption < < 1 implies that 1 + X/6 > a > X/5. Moreover, the assumption 

|(iy| > 1 implies additionally that a > 1. We can define the background stationary Markov 
process by 

Xt = RtY = sup {St - Ss- a{t - s)), t G R. 

— oo<s<t 

We have 

^ FtY. a — Xm 



«-A/[o,oo)^^(^^)' 

Unlike the previous example, here P{Xq = 0) = lim^joo Ee~^^" = a — X/S > 0. The local 
time L of X at has rate 

^ = i;'YiO+) =a- X/6. 

The assumptions on a imply that < 1 and hence we can construct the stationary process 
Q Qt = T^tL, where L(s, t] = L{s, t] — {t — s). We have 



P{Qo > x) = fie 



-e*x 



where 9* = tpY{9*) = X{a — 1) ^ — 6. Note that the latter is positive since a > 1 and 
1 + X/6 > a. 

4.3 Example 3: fluid queue driven by the local time of a risk-type process 

Let 

Yt = bt- St, t G M, 
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where 6 > 0, 5 is an a-stable subordinator, < a < 1, with 

and c is a positive constant. Thus, 5 is a (l/a)-self-similar process, i.e. {SKt,t € M) = 
{K}/°'St,t € M). We here have ESt = +oo for t > and S'f — > oo faster than Unearly, so 
Yt — > — oo, as i — > oo, a.s. Similarly, St — > — oo as i — > — cxd, a.s. So the stationary reflection 

of y 

Xt = ntY= sup {h{t-s)-{St-Ss)), t€M, 

— oo<s<t 

exists uniquely, due to Lemma [HI Physically, Xt is the content of a queue with linear input 
(arriving at rate h) and jump-type service represented by S. Alternatively, X is a so-called 
risk process in the theory of risk. We have 

iPy{0) = logSe-^(^*+i-^*) = he- cd"", 9>0. 

We refer to Section 12.11 and, specifically. Lemma [H for the distribution of Xq which is 
exponential with rate fi > where fj, satisfies = 0) 

^ = (c/6)V(i-°). 

The local time L of X at zero is such that t h-> L(0, t] is a.s. right-continuous (but not 
continuous) with rate /z. Assuming that /i < 1, or 

c < 6, 

we can further let L{s, t] = L{s, t] — {t — s) and let Q be defined by 

Qt = TZtL, t G M 

(see Lemma [TJ) Theorem [3] gives the distribution of Qq under Pq and under P. We have, 
PoiQo >x) = e-(^-^)^ P{Qo >x) = (0 ^ e-(^"^)^ x > 0. 

4.4 Example 4: fluid queue driven by the local time of a risk-type process 
with a Brownian component 

Take 

Yt = m + aBt- St, 

where S is the inverse local time of an independent Brownian motion. Assume o"^ > 0. We 
have 

ll;Y{e) = log ^e^(^l-^^) = m + - 2c0l/^ > 0, 

where c is a scaling parameter. Since limt_too Yt = — oo, a.s.. Lemma [8] allows us to construct 
Xt = TZtY . Here, Y is spectrally negative, and so, as shown in Lemma [T] 

P(Xo >x) = e"^^, a; > 0, 

where // > and 

V'y(/i) = 0. 
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Letting 

S = 1 + 5V"2c-2^ 

we find 

C ((5^ + 1)1/3 + (^7 _ 1)1/3 



Here, P{Xo = 0) = 0. As in Proposition [21 this /.i is the rate of the local time L of X. Note 
that, since Y has unbounded variation paths, the local time L is a.s. continuous. Assuming 
that < 1, which is equivalent to 

^y(l) =6+^Ct2-C>0, 

we construct Q as before: Qt = TZtL,t G M. From Theorem [3] we have that 

P{Qo >x) = ne-^^+^y-''^'', X > 0. 



5 Idle and busy periods 

In this section we study idle and busy periods of the fluid queue process {Qt : t G M} as 
defined in (jlSp . We work under the assumptions that Y is either spectrally negative or 
spectrally positive and that < < 1, where fi is the rate of L-see (flTl) . Under these 
assumptions, the process Q constructed above is stationary and the sets 

{t G M : Qt^> 0, Qt = 0} = {g{n) : n £ Z} 
{t£R: Qt-= 0, > 0} = {d[n) : n £ Z} 

are a.s. discrete, with elements denoted g{n), d{n), respectively. We need a convention for 
their enumeration, and here is the one we adopt. 
First, 

•••<(7(-l)<ff(0)<0<5(l)<5(2)<---, 

Second, 

g{n) < d{n) < g{n + 1), n € Z. 

Let A^i (resp. A''2) be the random measure (point process) that puts mass 1 to each of the 
points g{n) (resp. d{n)). Notice that the point processes Ni,N2 are jointly stationary. 

The intervals {g{n),d{n)) are called idle periods, while the intervals {d{n),g{n + 1)) are 
called busy periods. An observed idle period is, by definition, equal in distribution to an idle 
period, given that the idle period contains a fixed time t of observation. By stationarity, we 
may take the time of observation to be t = 0. In other words, 

observed idle period := ( {g{0),d{0)) \ Qq = O) = { {g{0),d{0)) \ g{0) < < d(0)). 

Here, = denotes equality in distribution under measure P. Similarly, 

observed busy period := ( {d{0),g{l)) | Qq > O) = ( {d{0),g{l)) \ d{0) < < ^(1)). (24) 

In this section we identify the distribution of observed idle (and busy) periods (see 
Propositions [5l [6|). In both cases we shall appeal to the result of Lemma [2] below a short 
proof of which is provided. Note that formula (j25p below and related facts are also proved 
in Kozlova and Salminen [101 Section 5] and Sirvio (Kozlova) |16j . 
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sample path of Q 




8(0) 8(1) 



Figure 1: Typical behaviour of Q and the background Markov process X when the underlying 
Levy process Y has unbounded variation paths. By convention, the origin of time is contained 
between g{0) and d{0). Note that excursions of X away from correspond to intervals over 
which Q decreases. 

Lemma 2. Let (0, P) be a probability space endowed with a P-preserving flow {9t, t € M) 
[see Appendix A]. Let Ni, N2 be jointly stationary simple random point processes (NioOi[B) = 
Ni{B + t), t £R, B e ^(M), i = l,2) with points {ti{n),n i = 1,2, such that 

■■■ < ti(-l) < ti(0) < < ti{l) <ti{2) <■■■ , 

and 

ti{n) < t2{n) < ti{n + 1), for all n e 1^. 

Let M be the random measure which is defined through its derivative with respect to the 
Lebesgue measure as 

M{dt)/dt = ^l(ti(n) < t < t2{n)). 

Assume that X\ has finite intensity. Let Pm be the Palm measure with respect to M . Then 

„.,<„,.,..<0)|,°^«r°'"°'l-f«b-'"""l_ p5) 

a — p 

Proof. It is easy to see that M is also stationary, i.e. MoOt{B) = M(B + 1). Let P/v- be the 
Palm measure with respect to Ni, i = 1,2 and let A be the intensity of A^i (which is-due to 
the law of large numbers-the same as the intensity of It follows easily from Campbell's 
formula that M has finite intensity: Am = A£^ArJt2(0) — ^i(O)] < 00. The Palm exchange 
formul^l between Pm and P/Vj yields 

rhil) 1-12(0) 

XmEm[Y] = XEn, / Yo6t M{dt) = XEn, / YoO^ dt, (26) 

Jii(O) JtiiO) 

^In [Tj p. 21 the formula is given and proved for point processes but the generalization for arbitrary 
jointly stationary random measures is straightforward. 
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sample path of Q 




Figure 2: Typical behaviour of Q and the background Levy process Y , in case that Y is 
spectrally negative with bounded variation paths. 

for any bounded random variable Y. Apply (p6]) with Y = e~"*2(o)+/3ti(o)^ Since ii(0) = 
sup{t < : Ni{{t}) = 1}, t2(0) = inf{t > ti{0) : N2{{t}) = 1}, and PiVi(ii(0) = < 
t2(0)) = 1, we have ti(0)o6'i = -t, t2{0)o9t = t2{0) - t, PN^-a.s. on {ti(0) < t < ^2(0)}. 
Therefore, Yo0t = e-<^t,io)^ia-i3)t ^ P^^-as. on {ii(0) < t < t2{0)}, and so 



(27) 



a — /3 

Arguing in a similar manner, through the exchange formula between M and we obtain 



(3 — a 



Setting /3 = 0, and then q = 0, in (j27p we obtain: 

l-E7Vi[e-"*2(°)]i 



a 



I - EN,[e-f'^(o^ 



/3 



On the other hand, with a = in (I28p . we have 



AMi5;A/[e-^*^( 



A 



1 - ENAe^''^^^] 
P 



(28) 

(29) 
(30) 

(31) 



The exchange formula between A^i and shows that the right hand sides of (jSOp and (I3ip 
are equal. Substituting these and (p9]) into ([27j) we obtain the result. □ 



5.1 Observed idle periods 

We are interested in the distribution of the idle period {g{0), d{0)), given that Qq = 0. 
The rationale used for this computation is as follows: We can always assume that Xq > 0, 
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sample path of Q 




Figure 3: Typical behaviour of Q and the background Levy process Y , in case that Y is 
spectrally positive with bounded variation paths. Here, only the case where the jump part of 
Y is compound Poisson is depicted. When Qt > and Xt = 0, we see that Qt increases at 
rate |(iy| — 1, where dy < —1 is the drift ofY. 



since this is an event with probabiHty 1. If Qo = then Qt wiU remain at least until X 
hits 0, since Q cannot increase unless there is an accumulation of local time L, and this 
can happen only when X is 0. Recall that the first hitting time of by X is denoted by 
D = inf{i > : Xt = 0}. If Qq = 0, the first time that Q becomes positive has been denoted 
by d{0). Our claim is: 

Lemma 3. Given that Qq = the ending time of the idle period is a.s. equal to D, i.e., 

P((i(0) = D\Qo = 0) = l. 

Proof. From the argument above we have that D < d{0) a.s. on {Qq = 0}. Suppose that 
there is JIq C with P{^o) > such that Qq = and D < d{0) a.s. on Qq. If Qq = and 
D < d{0) then 

Qt = sup {L{u, t]-{t- u)} = 0, for all t € {D, d{0)). 

D<u<t 

This implies that 

L{u, t] <t - u, for all D < u <t < d{0), 

which means that, for uj € Qq, L{uj,-) is absolutely continuous on some right neighbourhood 
of D. If Y is spectrally negative or if Y is spectrally positive but not of bounded variation, 
then L is a.s. singular on any right neighbourhood of D, and we obtain a contradiction. If 
Y is spectrally positive with bounded variation paths then L is absolutely continuous and 
is given by (fT3|) . In this case, Q increases at rate |(iy| — 1 > whenever it is positive and 
this shows immediately that here, too, D = d{0) a.s. on {Qq = 0}. □ 

Remark 1. The result in Lemma[3]can alternatively be expressed by saying that the process 
{Lt — t : t > 0) is under Pq initially increasing. In [13] Proposition 6.3 this is proved in 
the case X is a reflecting Brownian motion, and the proof therein could have been modified 
to cover the present case. However, we found it motivated to give the above proof which 
highlights other aspects than the proof in [13] . 
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Using Lemma [2l we shall reduce the problem to that of finding the distribution of 
D = inf{t > : Xt = 0} given that Qo = 0. Let A''i (resp. N2) be the point process with 
points {gin} (resp. {d{n)}). Then M{dt)/dt = ±{Qt = 0), and so 

Pm = Pi- 1 Qo = 0). 
Formula (j25p . together with Lemma O then gives 

E^^-am+MO) I Q, = o] = aE[e-^''\Qo = 0]-PE[e-nQo = 0] _ ^32) 

a — P 

To compute the distribution of D given Qo = we need the following two lemmas. 
Lemma 4. Let 

G := sup{t < : Xt = 0}. 

Then it holds 

{Qo = 0} = {Qg + G<0}. 
Proof. Since Xt > for all t G (G, D), we have 

L{s,t]=0, G<s<t<D. 

Recall that 

Qt = ns,tL{Qs) = sup L{u,t]y {Qs + L{s,t]) 

s<u<t 

So, if G < s <t < D, we have L{s,t] = L{s,t] — (t — s) = —{t — s), i.e. 

Qt = iQs-it-s))+, G<s<t<D. 
If we assume that Qo = 0, we have G < g and so 

= Qg = {QG-{g-G))+, 

which implies that Qg + G = g < 0. □ 

Lemma 5. (i) Conditional on Xq, the random variables Qg,G,D are independent (under 
P)- 

(ii) For allx>0,t> 0, 



P{Qg >t) = P,{Qg >t)= PoiQo >t) = e- 



rt 



where 6* = V'y(l) ^/ ^ spectrally negative or is equal to the unique positive solution of 

6* = iPy{G*) ifY is spectrally positive. 

(Hi) Qg is independent of{G,D) (under P). 

Proof, (i) The independence follows from the strong Markov property at G (at which 
Xg = 0) and Markov property at 0. Indeed, first observe that G is a stopping time with 
respect to the filtration {^t := it(X_s,0 < s < t), t > 0}. Second, Qg = 'R-gL = 
suPs<G ^{s, G] = supg<(^ [L{s, G] — {G — s)) and so Qg^{G < t) is measurable with respect 
to := a{Xs,s > t) for all t. This proves independence between Qg and G. Third, D 
is measurable with respect to = a{Xs,s > 0). So, conditionally on Xq, the random 
variable D is independent of the pair {Qg,G). (ii) The distribution statement about Qg 
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follows from the strong Markov property at G. Let, as usual, = £ > 

0): An {-G <t}£ ^_ J. Since Qg = QqoOg, 

P{Qg >t) = PiQooOG >t) = EPiQ^oOG > t I ^_g) = EPxaiQo > t) = Po(Qo > = e"^** 

where the latter follows from Theorem [3l (hi) This is immediate from (i) and (ii). □ 
Proposition 5 (distribution of observed idle period). Fix a,(3 >0, (3. 

(i) When Y is spectrally negative we have 

$y(0) V'y(l) f a $y(a) - 1 /3 $r(/3) - 1 



1-$y(0) a-P \a-ijY{l) ^>y(a) /3-Vy(l) ^y{P) 

(ii) When Y is spectrally positive we have 

-„,(o)+/35(o) in - nl - ^* f a-<^Y{a) _ /3-<^y{(3) 

^ J 1-V'y(0+) a-/3Va-^* 

w/iere 6* > is defined by ipyiO*) = 

Proof. From Lemma [5] we have that D,Qg are conditionally independent given Xq and G. 
Hence 

£[e-^^l(QG + G < 0) I Xo, G] = i?[e-^^ | Xq, G] P{Qg < -G \ Xq, G). 

= E[e~'^\Xo,G] (1-e^*^) 

where the second line was obtained from the facts (all consequences of Lemma [5]) that (i) 
D, G are conditionally independent given Xq, (ii) Qg, G are also conditionally independent 
given Xq, and (iii) Qg is independent of Xq and exponentially distributed with parameter 

0* = /^^(^)' ^ spectrally negative ^^^^ 
1 ■i/'y(0*), if Y is spectrally positive 

Taking expectations we get 

E[e-'^1{QG + G<0)]= i?[e-^^ - e"^^+^*^], 

and, using Lemma HI 



^[e-*^^ I Qo = 0] 



i^(Qo = 0) 



We now use a version of LemmalU formulated for excursions of general stationary processes; 
Pitman [141 Corollary p. 298; references therein]. The random measures Ni,N2 correspond 
to the beginnings and ends of excursions of the stationary process {Xt,t € M), and, since 
there is never an interval of time over which X is zero, the random measure M coincides with 
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the Lebesgue measure, while Pm = P- Applying Pitman's result-an analogue of formula 
([25|) -gives 

a — (3 

The joint Laplace transform of D,G is thus expressible in terms of the Laplace transform 
of D. Combining the last two displays we obtain: 



Using this in (j32p results in 

1 



^[g-ad(0)+/3c;(0) | = Q] 



P(Qo = 0) a -(3 

X { (i? [e-^*^] - E[e~-^]) - ^ (i? [e-^*^] - E[e-^^]) } (35) 

So far, the arguments are general and hold for both spectrally negative and positive Levy 
processes y, as long as 9* is taken as in (ISSp . Substituting next the expression for the 
Laplace transform of D from ([8]), (fT6l) for the spectrally negative, respectively positive, 
case, we obtain the result. □ 



5.2 Observed busy periods 

In this section we follow ideas in [TB] . We are interested in the distribution of the observed 
busy period, as defined in (|24p . On the conditioning event {Qq > 0}, we have, by our 
enumeration convention, 

5(0) < (i(0) < < 5(1), P-a.s. 

Using Lemma [2] with A''i (resp. A''2) the point process with points {d{n)} (resp. {g{n}), we 
have 

i;[.-«w+«o) I Q„ > 01 = °i^b-°-""IQ.>o]-gE[^-^""iqo>o|^ „ ^ „ (3,) 

a — (3 

Recall the evolution equation for Q: 

Qt = Qs + L{s, t]-{t-s)- inf {Qs + L{s, u] - {u - s)}. (37) 

s<w<i 

Let s = and assume Qo > 0- Since Xq > 0, P-a.s., we have L(0, t] = for all < t < L> = 
inf{r > : Xr = 0}, and so, 

Qt = Qo - t - inf {Qq - u} = Qo - t, a.s. on {Qq > 0, t < D}, 

0<u<t 



which implies that 



5(1) = I {0 < Qo < D}, ^3g^ 

\g{l)oeD, a.s. on {Qo> D}. 
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Now, if Qo > D, we have Qd- = Qo — D, so from (i37jl . Q evolves as 

QD+t = Qo- D + L[D,D + t]-t, t>0, 
ans as long as Qo+t > 0. This implies that, a.s. on {Qq > D}, 

giiyOo -D = inf{i > : Qq - D + L[D, D + t]-t = 0}. 
Therefore (l38l) becomes 



\D + mi{t>Q: QQ-D + L[D,D + t]-t = Q}, a.s. on {Qo > 

Consider now the inverse local time process, with the origin of time placed at D, i.e. 

L^i^ := inf{t > : L[D, D + t\> x}, x > 0. 

By the strong Markov property for X at the stopping time D we have that the P-distribution 
of {L~^^^, X > 0) is the same as the Po-distribution of {L~^, a; > 0), which has been identified 
in Propositions [2] and m Thus, {L~j^^, x > 0) is a (proper) subordinator. Consider next the 
spectrally negative Levy process 

Aj, := X - L^l^, X > 0. 

Notice that P(Ao = 1). The Laplace exponent of A is the function tjj\ of \22\\ . Define the 
hitting time of level —a by A 

o-(A; a) := inf{x > : A^; < -a}, a > 0, 

Formula (I6ip gives us the Laplace transform of o"(A; a) in terms of the scale functions of A, 
defined in (j56p and (j57|) . Combining them, we obtain 



As can be easily seen from Lemma [H for any a > 0, 

inf{t > : t - L[D, D + t]>a} = inf{x > : A^. < -a} + a = a(K; a) + a. 
Using this in (I39p . we obtain 



5(1) 



Qo, a-s. on {0 < Qo < D}, 

Qo + (j(A;Qo - D), a.s. on {Qo > D}. 



It is useful to keep in mind that A is independent of Qo ~ -C) by the strong Markov property 
of X at D. We are now ready to compute the Laplace transform appearing on the right 
hand side of (l36l): 



E[e~'^9{i) . > 0] = ^[e-°f ; < Qo < P*] + E[e~'^9W -Q^y D] 

= P[e~"^«; < Qo < + P[e-°(^o+'^(^;^"-^); Qo > D] 

= P[e-°^°;Qo > 0] - E[e~''^° (1 - e-"'"(^'^°-^)); Qo > D] (41) 
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Recall that P{Qq > x) = fie ^*^, and so 

£;[e-"«";Qo>0]=M^^. (42) 

To compute the second and the third terms we need some elementary properties of expo- 
nentially distributed random variables which we state without proof. 

Lemma 6. Let T he an exponentially distributed r.v. with parameter A and {X, Y), X >0, 
Y >0, be a two-dimensional r.v. independent ofT. Then X andT — X — Y are independent 
given T > X + Y . Moreover, 

S[e-"(^-^-^)|T>X + y] ^ 



a + A' 

£;[e-°^; T>X + Y]= ^[e-("+^)^-^^]. 

Use ([37|) once more with s = G = sup{t < : = 0}, and t = 0, taking into account 
the fact that L is not supported on (G, 0), to obtain 

Qo = Qg + G, a.s. on {Qo > 0}. 

Since Qg is exponentially distributed with parameter 9* and independent of (G, D) (from 
Lemma [5]), we have, applying Lemma [6j the following result: 

Lemma 7. Given Qq > D, the r.v. 's Qq — D and D are independent. Moreover, 

i?[e-"W»-^)|Qo>l?] = -^, 

a + 0* 

E[e~''^; Qo>D]= iJ[e-(-+e*)^+e*G]_ 
Using Lemma [71 we write the last term of (|4ip as follows: 

= P{Qo > D) ^[e""^ e-°(Q«-^) (1 - g-^'^^^^Oo-^)) \Qo>D] 

= P{Qo > D) ^[e""^ \Qo>D] S[e-"Wo-^) (1 - e-"-(^'«o-D)) | Qo > ^] 

where, in the last term, we introduced a random variable V, exponentially distributed with 
parameter 6*, independent of everything else (due to the fact that Qo — D, conditionally on 
being positive, is exponential with parameter 9* , independent of A). The first term of ()43p 
can be computed as in ([M]) . We have, for all a, (3 > 0, a P, 



^ ^ a-/?Uy(a) 



$y(/3) 

and for the spectrally positive case, for all a, /3 > 0, a ^ (3, 

E[e-^^^<^]=M^+)^^^^^^^^^ 

a — p 
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Note that taking account of the definition ()22l) of ?/^a for both the spectrally negative and 
positive cases, and the fact that iIj\{9*) = 0, one sees that generically for both cases, for all 
a > 0, 

= /^(c, _ ^^(a + e*)). (44) 
a 

For the second term of (j43]) we have, using (j40]) . 



a + 

'0A(a + ^'*) a 



a + d* ipxia + 0*) - a \ a + 6* $A(a) 

~ a + 9* <I>A(a) i/;A(a + 6**) - a 
Multiplying ((44jl and (j45]) we obtain the following expression for (l43]l : 

a + y* 4*A(a) 



This, together with and (jH]), yields: 

£;[e-"s(i) I Qo > 0] 



^>A(a) 

Using p6p . we finally come to rest at the following main result. 
Proposition 6 (distribution of observed busy period). For a,f3>0, a ^ P, 

^ ^ a-/?UA(a) <^a{P)) 

where ^\ is the right inverse of ip^ which is given in \2S^ and ( cf. Theorem 0) 

(i) 6* > is defined by ?/'y(l) the case that Y is spectrally negative, 

(ii) 9* > is defined by ^y(0*) = 9* in the case that Y is spectrally positive. 



(45) 



5.3 Typical idle and busy periods 

We now consider the problem of identifying the distribution of a typical idle and a typical 
busy period of Q. We place the origin of time at the beginning of such a period, by 
considering the appropriate Palm probability. Let A'^ (resp. N^) be the point process with 
points {g{n)} (resp. {d{n)}), the beginnings of idle (resp. of busy) periods, and let Pg (resp. 
Pd) be the Palm probability with respect to Ng (resp. Pd)- 

X 

8(0) diO) 8(1) 
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Using ([29]) we have 

i^Q^[e-"'^(°)|Qo = 0] = 1 - Eg[e-'"^^% 
A 

where ^ = P{Qo > 0) and A is the common rate of A'^ and N^. The right side is precisely 
what we need. Everything in the left side is known (see Prop. [5]) except the rate A. Consider 
first the case when Y is spectrally negative. Using Prop. [5^i) with /3 = we get 

l-^^HO^ a ^y{a) - 1 _ ,(0) 

Taking limits as a — > oo-and since <I>y(a) oo-we find the value of A and so the Laplace 
transform Eg[e~""^^^^ of the typical idle period. The result is in Proposition [7Ki) below. 

We repeat the procedure for the spectrally positive case and, using Proposition [5^ii), we 
obtain: 

Note that limg_+oo ^y(^)/^ = oo if y is of unbounded variation, and so liniQ^oo = 
0. Thus, we can find A and £'g[e~"'^^''^]-see Proposition ©(ii) below. 

But if Y is spectrally positive (with non-monotone paths) and of bounded variation then 
L is absolutely continuous and has a drift dy-see (I13p . We can easily see, e.g. from (I12p . 
that 

poo 

^I^Yie) = iogi?[e-^(^^-^o)] = \dY\e - / (1 - e-'ymdy), 

Jo 

and, since f^{yAl)Il{dy) < oo, we obtain limg^oo '0y(^)/^ = \dY\- So lima^oo ^y(a)/a = 
l/|dy|. Again, we can find A and £'g[e~"'^('')]-see Proposition ©(iii) below. 

Proposition 7 (distribution of typical idle period). Fix a > 0. Let Pg be Palm probability 
with respect to the beginnings of idle periods of Q. 

(i) When Y is spectrally negative we have 

A = a>y(0)^y(l), Eg[e--^^^)] = 1 - 

$y(a) a - V'y(l) 

(ii) When Y is spectrally positive and yll{dy) = oo, we have 

A = ^Wo)r , Eg[e~-'(^^] = 1 - r 

where 9* > satisfies 6* = ipYi^*)- 

(iii) When Y is spectrally positive and Jq yll{dy) < oo, we have 



where 9* > satisfies 9* = ipYi9*), and dY is the drift defined in (fT2|) - (fT3|) . 

Remark 2. By Assumption [A2]-see Section [3l-we have (iy < — 1 and so the constant above 
is positive. 
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In the same vein, we obtain the Laplace transform £"^[6""^'-^-*] of a typical busy period. 
Notice that, under P^, we have d{0) = 0, and so the first busy period to the right of the 
origin of time is the interval {d{0), g{l)). We have, 

^aii;[e-"s(i)|Qo > 0] = 1 - Ed[e-"3ii)^. 



Using Proposition E] with /? = 0, we have 

l-Ed[e 



A <I>A(a) 

From the expression (p2|) for ij^AiO) we find limg^oo V'a(^)/^ = 1- So, limQ,^oo ct/^A(cK) = 1- 
Therefore: 

Proposition 8 (distribution of typical busy period). Fix a > 0. Let he Palm probability 
with respect to the beginnings of busy periods of Q. Let be defined as in (I22p . and let <I>a 
be its right inverse function. Then 

"9(1)1 - " 



A = fie*, Ed[e-''3(^^ = 1 



$A(a)' 



where fx = $y(0), 6* = ^/;y(l), ifY is spectrally negative; and /x = ^/;y(0), 9* > {) is defined 
through 6* = ■0y(^*), ifY is spectrally positive. 

Corollary 1. The mean duration of a typical idle period is (1 — fj,)/X, while the mean dura- 
tion of a typical busy period is fi/X, where fi is given by (fTTI) and A is given in Propositions 

Remark 3. By the relation between P and the Palm probability P^ it follows that the 
typical idle (respectively busy) periods are stochastically smaller than the observed idle 
(respectively busy) periods, see [T]. In particular, the means of the former are shorter than 
the means of the latter. (This is usually referred to as the "inspection paradox".) This gives 
us several inequalities between different quantities associated with the Levy process Y. To 
give an example, we compare the mean durations of idle periods in the spectrally negative 
case. From Proposition [5] we have 

L J l-fy(O) Uy(«)(«-^y(l))y 

_ $y(0)^y(l) 



l-$y(0) 

It follows that 



F{a). 



M''!") I * ^ 0] . -MM no). 

Since d{0) and —g{0) are identical in law, the mean duration of the observed idle period is 

£[d(0)-9(0)|Q„ = 0]=-2M^f'(0). 

Notice that, using the function F, we may write from Proposition [7] 

Eg[e-^m] = l-aF{a). 



26 



and 

Eg[dm=F{0). 

Now, by the "inspection paradox" , 

Eg[dm <E[diO)-giO) \Qo = 0], 
which, after some manipulations, is equivalent with 

(1 - ^Y{0)f < 2 ($y(0)2 - «.y(0) + ^^^^(l)) • 

Example 5 (continuation of Example 1): Consider Yt = aBt — fit, and assume < 
fi < 1. Here the rate of beginnings of idle (or busy) periods is 

The mean duration of a typical idle period of Q is 

2/x' 

while the mean duration of a typical busy period of Q is 



2(1 -/i)- 

To find, e.g., the distribution of a typical busy period, we use Proposition [HI We have, see 

i'kiq) = q- '^riq), 

where is the inverse function of V'yj i-e. 



and <I>A is the inverse function of tpA, i.e. 

^ , , (1 - ^) + 2a^a + J{1 - /i)2 + Aa^a 
^A{a) = , 

and so the Laplace transform of the typical busy period is 



(1 - ^) + 2a^a + V(l-^)2 + 4cj2a' 

A On Skorokhod reflection, fluid queues, and stationarity 

In this section we review some facts about the Skorokhod reflection of a process with 
stationary-ergodic increments. We carefully define the system, give conditions for its stabil- 
ity, and recall some distributional relations based on Palm calculus, see [I], [8]. Although 
the setup is much more general than the one used in later sections for concrete calculations. 
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it is nevertheless interesting to isolate those properties that are not based on specific dis- 
tributional assumptions (such as Markovian property or independent increments) but are 
consequences the more general stationary framework. 

Let {Q, P) be a probability space together with a P-preserving flow {6t, t G R). That 
is, for each t € M, 0* : O ^ O is measurable with measurable inverse, is the identity 
function, OtoOs = 6s+t, for all s,t eM, and PiOtA) = P{A) for alH G M, ^ e ^. Consider a 
process W = (Wj,t G M) with stationary increments, i.e. (Wj — \Ns)o0^ = Wi_|_„ — \I\/s+u for 
all s,t,u G M. We let E denote expectation with respect to P. 

Following [8], we define the "Skorokhod Dynamical System" (abbreviated SDS henceforth) 
driven by \N as a 2-parameter stochastic flow: 

7^,,^W(x) := [x + - W,] - inf [{x + Ws- Wu) A 0] (46) 

s<u<t 

:= sup (Wt - W„) V (x + Wt - W,) X > 0, s<t. (47) 

s<u<t 

Thus, for each s < t, we have a random element TZg^t^ taking values in the space of C(M+) 
of continuous functions from into itself. The family (T^-^^jW, — oo < s < i < oo) is a 
stochastic flow because the following composition rule (semigroup property) holds for each 
w G r2: 

7^,,tW = 7^„,^Wo7^,,„W, s<u<t, 
nt,t\N{x)=x, teR,x>0, 

It is a stationary stochastic flow because, for each x G M-|-, we have: 

7^,,tW(x)oe„ = 7e,+„,t+„W(x), -oo<s<t<cx), mGE. 

We say that the process Z = jZ^, t G M} constitutes a stationary solution of the SDS driven 
by W if Z is W-measurable and if 

Zt = 7^,,tW(z,), s < t, 

Existence and uniqueness is guaranteed under some assumptions: 

Lemma 8. Assume that sup_oo<^<oWs < oo, and lim ^_,^ \Nt < oo, P-a.s. Then there is 
a unique stationary solution to the Skorokhod dynamical system driven by W. This is given 
by 

Zt= sup {\Nt-\Nu) =:nt\N. (48) 

—oo<u<t 

Quite often, in addition to stationarity of the flow, we also also assume ergodicity, namely 
that each A ^ that is invariant under 9t for all t, has P{A) equal to or 1. Owing to 
Birkhhoff's individual ergodic theorem Lemma [8] immediately yields: 

Corollary 2. Under the ergodicity assumption, and if E\Ni < 0, then there is a unique 
stationary solution Z to the Skorokhod dynamical system driven by W. The process Z is 
given by (148 p and Z is an ergodic process. 
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For the purposes of this paper, assume that W is of the form 

\Nt-\N, = A{s,t]-f3{t-s), s<t, 
where ^4 is a locahy finite stationary random measure, and 

0< f3 <a:= EA{0,1). 
Let Pa be the Palm probabihty, see [I], [6], with respect to A: 



(49) 



(50) 



1 



PAiC) = -E 



a 



(0,1] 



±c°OtA{dt) 



The following is a consequence of Theorem 3 of [S]: 

Lemma 9 (distributional Little's law). Let Z he the unique stationary solution to the SDS 
driven by W of the form ([19|) . Assume that ([50j) holds. Then, for any function ip : [0, oo) — > 
M, which is continuous on (0, oo), we have (Theorem 3 of J3i) 



1 - ^(0) + ^EA^PiZo). 



(51) 



In particular, 



It should be noted that the decomposition (j49|) of W is not unique; nevertheless, ([ST 
holds, regardless of which decomposition of W we choose. 



P{Zo >x) = -j^Pa{Zo > x), x> 0, P(Zo > 0) = ^. 



B Exit times for spectrally negative Levy processes 

In this section we consider a spectrally negative Levy process and some facts regarding the 
first time the process exits an unbounded interval. Let Y = (Yt,t G M) be a spectrally 
negative Levy process and Levy measure 11. In other words, let B he a standard Brownian 
motion, r] an independent Poisson random measure on M x M_ such that 

Er]{dt,dy) = dtU{dy), U{0} = 0, / {y^ A l)U{dy) < oo, (52) 

let a G M, (7 > 0, and define, for —oo < s < t < oo, 

Y{s,t]=a{t-s) + a{Bt-Bs) 

+ / / yr]{du,dy)+ / y [r]{du,dy) - dull{dy)]. (53) 

J(s,t] J(-oo,-l] J{s,t] J{-lfl) 

Notice that we have thus defined only the increments of Y; the exact value of Yq is unim- 
portant; we may, arbitrarily, set 

Yo = 0. 
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(The reason that increments are more fundamental than the process itself is amply explained 
in Tsirelson [IT].) If we set 



y.:=r''*]' tern, 

\-Y{t,0], t<0 

we have 

Y{s, t]=Yt- Ys, -oo < s < t < oo. 

(In case that f^-^ \y\Il{dy) < oo, and o" = 0, the process Y has bounded variation paths and 
can also represented as 



Y{s,t]=dYit-s)+ / yr,{du,dy), (54) 

J{s,t] J(-oo,0] 

for some constant dy known as the drift of y.) 

To be more precise, especially for the construction of the stationary versions of processes 
in this paper, we introduce shifts. Assume that {B, rj) is defined on a probability space 
(r2,^,P) taken, without loss of generality, to be the canonical space Q. = C(M) x 7V(M^), 
where C(M) are the continuous functions on M, and A/YM^) are the integer-valued measures 
M?. Let P be the product measure on the Borel setqfl of C(R) x A/'(M^) that makes B a 
standard Brownian motion and r] a Poisson random measure with mean measure as in (I52p . 
and to each lo = (</?,//) € C(M) x 7V(M^), let B{ip,ij,) = (p, r]{ip,iJ,) = ji. Consider also the 
natural shift {Ot,t G M) on Vl defined by 

et{^,^l){s,A) = {^{t + s),^l{A + s)), sgm, A^m{M^), 

where A + s := {{t + s,y) G : {t,y) € A}. By construction, Y has cadlag paths, and, 
under P, it has stationary (and independent) increments. Henceforth we shall denote by P^ 
the conditional probability of P given Yq = x and expectation with respect to it. All of 
the following facts are standard results which can be found, for example, in [2l [11] See also 
|12j for a review which is more convenient for the setting at hand. 

Let : M ^ C denote the characteristic exponent of Y: 

and let tjjY '■ [0, oo) i-^ M denote the Laplace exponent of Y: 

V'y(/3) =logi?[e^^i], p>0. 

It is well known that ipy is infinitely differentiable, strictly convex, ^(0) = 0, lim^^oo (/?) = 
oo, and 

^^.(0+) = EYi = E{Yt+i -Yt) eRU {-oo}. 

For each (7 > let 

^Y{q) = sup{(3>0:i>Y{P)=q}. (55) 

Since Y drifts to infinity, oscillates, drifts to minus infinity accordingly as V'y(0+) > 0, 
iP'y{0+) = and iP'y{0+) < 0, it follows that $y(0) > if and only if Vy(0) < and 
otherwise <^y(0) = 0. It is also easy to see that ^'y((7) > for all q > 0. 



The space J\f{M.'^) is endowed with the topology of weak convergence; see, e.g., [7]. 
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Define also the scale functions W^'^\x), Z^'^\x) via their Laplace transforms 



oo 

e~ 







oo 



e 



-I3x 



defined for all p > ^riq)- The functions $y, 1^(9), appear in the expressions for the 
Laplace transform of the first passage times 

r+ := inf{t >0:Yt> x}, (58) 
tZ^ = inf{t > : < -x}, (59) 

as follows (cf. Chapter Theorem 8.1, p214 of Kyprianou |llj): 
Lemma 10. For all q > 0, x > 0, 

= e-*^('')^ (60) 
^[e""--.] = Z^i^x) - -^^^^'^(a^)- (61) 
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